Experimental time-resolved photoemission and ab initio study 
of lifetimes of excited electrons in Mo and Rh 



O 

o 



(N 






c^ 



O 

(N 
> 

\o 

O 

o 
\o 
o 



o 



X 



A. Monnich, J. Lange, M. Bauer, and M. Aeschlimann 
Department of Physic, University of Kaiserslautern, 67663 Kaiserslautern, Germany 

I. A. NechaeviB V. P. Zhukov\ P. M. Echenique^'^ and E. V. Chulkoyi'^ 

^ Donostia International Physics Center (DIPC), P. de Manuel Lardizahal, 

4, 20018, San Sebastian, Basque Country, Spain 

^ Departamento de Fisica de Materiales, Faeultad de Ciencias Quimicas, 

UPV/EHU and Centro Mixto CSIC-UPV/EHU, Apdo. 1072, 20080 San Sebastian, Basque Country, Spam 

(Dated: February 6, 2008) 

We have studied the relaxation dynamics of optically excited electrons in molybdenum and 
rhodium by means of time resolved two-photon photoemission spectroscopy (TR-2PPE) and ab ini- 
tio electron self-energy calculations performed within the GW and GW -\- T approximations. Both 
theoretical approaches reproduce qualitatively the experimentally observed trends and differences 
in the lifetimes of excited electrons in molybdenum and rhodium. For excitation energies exceeding 
the Fermi energy by more than 1 eV, the GW -\- T theory yields lifetimes in quantitative agreement 
with the experimental results. As one of the relevant mechanisms causing different excited state 
lifetime in Mo and Rh we identify the occupation of the 4d bands. An increasing occupation of the 
Ad bands results in an efficient decrease of the lifetime even for rather small excitation energies of a 
few 100 meV. 

PACS numbers: 71.15.-m, 78.47.-t-p, 79.60.-i 



I. INTRODUCTION 

The dynamics of excited electrons play an impor- 
tant role for many physical and chemical phenomena. 
Detailed knowledge of relaxation times of excited elec- 
trons is substantial for the interpretation of photoemis- 
sion spectra, surface chemistry data and low-energy elec- 
tron diffraction data. For instance, the relaxation time 
determines the inelastic electron mean free path, and 
hence the escape depth length of the applied surface 
analytical technique.^ Particular in the optical excita- 
tion regime {E < -Evac) experimental data are very 
rare due to the limited experimental access to this en- 
ergy regime. Therefore, the knowledge of hot-electron 
lifetimes is still far from being complete. Among ex- 
perimental methods developed for such studies one of 
the most suitable is the technique of time-resolved two- 
photon photoemission spectroscopy (TR-2PPE), which 
enables lifetime measurements directly in time domain 
at a temporal resolution of a few femtoseconds^. In the 
last decade several TR-2PPE experiments have been per- 
formed for simple metals'^ and noble metals^'^iSiLS fer- 
romagnetic 3(i metals (Fe, Co, and Ni))24i& and high-Tc 
superconductoraii. Of the paramagnetic transition met- 
als, the electron lifetimes have been measured only for 
the 5d-metal tantalumJ^. However, so far, no system- 
atic experimental studies within the Ad and 5d-transition 
metal series have been undertaken. 

Theoretical calculations of excited electron lifetimes 
have been performed in the past within the GW ap- 
proximation (GWA) for electron self-energy for bulk 
nobleiii^ii^ and Ad transition metalsi^^iSiiiiSiiS as well 
as for hd transition metal Taj^S In contrast to noble met- 



als, which show qualitatively similar band structure and 
density of states (DOS),^" electronic structure of Ad met- 
als varies strongly on moving from the beginning of the 
4(i-series to the endi^ The calculations by Zhukov et ali^ 
and Bacelar et ali^ showed that the evaluated lifetimes 
also strongly vary along the 4(i-series following trends 
in electronic structure. The extension, GW + T, of the 
GW approximation by inclusion of multiple electron-hole 
scattering within a T-matrix approximationSLSSiSiSiiS^ 
results in the decrease of the GW lifetime value that 
brings theory and experiment to better agreement ^^jSaS^ 



In this paper we perform TR-2PPE measurements of 
excited electron lifetimes in Mo and Rh and compare 
these results with the GW and GW + T calculated life- 
times. The choice of these metals is motivated by their 
very distinct electronic structure. The Fermi level (Ep) 
of Mo, a typical Ad transition metal with bcc crystalline 
structure, is located in a deep minimum of a wide and 
half-filled d-band (see Fig. 4)iS£ Rhodium in contrast, a 
Ad metal with fee crystalline structure, exhibits a very 
narrow and almost-filled d-band^, where the Fermi level 
is situated in the vicinity of a sharp peak of density of 
states near the edge of the d-band (see Fig. 7). Hence, 
in comparison with molybdenum, d states of Rh are spa- 
tially more localized and, consequently, multiple scatter- 
ing contributions to the electron lifetime should be en- 
hanced. Such considerable differences in the band struc- 
ture between Mo and Rh are expected to cause strong 
distinctions in the excited state lifetimes. 




FIG. 1: 2PPE cross-correlation trace at an intermediate en- 
ergy E — Ef — 1.3 eV; inset, principle of a TR-2PPE pump 
and probe experiment. 



II. EXPERIMENT 

With the TR-2PPE technique the decay dynamics of 
an electron excitation (intermediate state) located be- 
tween the Fermi edge Ep and the vacuum level E'vac can 
be studied. The intermediate state is populated by a first 
(pump) laser pulse. Due to electron-electron scattering 
processes the excited electrons will then decay to lower 
unoccupied states. The rate of inelastic scattering of the 
electron defines the inelastic lifetime of the intermediate 
state. After a controlled delay t a second (probe) laser 
pulse transfers the remaining population to the detected 
state above the vacuum level. The raw data (see Fig.^ of 
a TR-2PPE measurement exhibits a 2-pulse correlation 
trace at a fixed kinetic energy, where the shape of the 
trace contains information about the pulse width of the 
pump and the probe laser pulse as well as the population 
decay time (and hence inelastic lifetime) of the probed 
intermediate state. In order to extract the intermedi- 
ate state lifetime from the crosscorrelation trace, a rate 
equation model of a three level system is solved within a 
matlab-based fitting routine which uses a least-square-fit 
algorithm. Such a rate-equation model corresponds to 
solutions of the Liouville-von Neumann equations within 
the density matrix in the limit of rapid dephasing. It has 
been shown before to be a equivalent (and efhcient) ap- 
proach for the description of electronic band excitations 
as investigated in the present studies^^-. 

The laser pulses for the TR-2PPE experiments were 
delivered by a mode-locked Ti:sapphire laser, pumped by 
about 7 W from a cw NdiYVOa diode laser. The available 
energy range for 2PPE emission has been increased by 
doubling the frequency of the linearly polarized pulses in 
a 0.2 mm thick beta barium borate (BBO) crystal. 

To minimize coherent effects in the 2PPE signal the 
polarization of pump and probe pulse were set to s and p 
polarized (cross polarized), respectively. The coUinearly 



recombined beams were focused at 45° incidence with re- 
spect to the surface normal of the sample. With these 
UV pulses {hi' — 3eV) of up to 0.75 nJ, at an repetition 
rate of about 80MHz less than one electron per pulse was 
emitted from the sample surface. In this means we avoid 
space-charge effects which can significantly distort the 
2PPE signal. Furthermore it guarantees that the excita- 
tion density of the intermediate excited state is kept at a 
low level so that the experiment probes the relaxation of 
individual excited electronic states rather than the collec- 
tive behavior (thermalization dynamics) of a transiently 
heated nonequilibrium distribution. A cylindrical sec- 
tor energy analyzer (CSA) — positioned in the direction 
of the surface normal (/cy = 0) — was used to detect the 
electrons with an energy resolution of 100 meV. 

For the non-time-resolved 2PPE measurements, 
we used a cw argon-ion-laser pumped mode-locked 
Ti:sapphire laser delivering transform-limited pulses of 
about 120 fs, 22.5nJ/pulse at a repetition rate of 80 MHz. 
The advantage of this narrow bandwidth laser system is 
its tunability in the laser wavelength within the amplifi- 
cation bandwidth of about 100 nm which is required to 
distinguish between spectral features in the 2PPE spec- 
tra, arising from the occupied and unoccupied density 
of states of the electronic structure (see below). •^^ The 
wavelength of the laser pulses can be varied from 740 to 
840 nm which corresponds to a tuning range of 3.35 to 
2.95 eV in the photon energy of the frequency-doubled 
laser pulses. 

Metallic Mo and Rh sheets of approximately 0.1-mm 
thickness were used for our experiments. We chose poly- 
crystalline samples to guarantee an averaging over all 
possible directions in the k space of the Brillouin zone 
and to avoid any influence of surface resonances on the 
measured lifetime which can exceed the bulk lifetime val- 
ues by a factor of 10 and moreiSl The purity of the sam- 
ple was specified to be higher than 99.9%. The metallic 
sheets were cleaned by several sputter and heating cy- 
cles. After the cleaning we found no indication of sig- 
nificant impurities, such as C or O at the surface. Next 
to the sheet samples further TR-2PPE experiments were 
performed at ultrathin polycrystalline films of molybde- 
num of varying thickness between 20 nm and 100 nm. 
The films were prepared in-situ by means of electron- 
beam evaporation on a dielectric MICA substrate from 
a molybdenum rod (purity, 99.9-|-%) at a pressure of 
< 10~^ mbar. The film thickness was controlled using 
a calibrated flux monitor. 

The cleaned metal surface was dosed with a small 
amount of Cs (< 0.1 ML) to lower the surface work func- 
tion. This procedure enables to increase the accessible 
energy range in the 2PPE experiment by about 1 eV so 
that excitation levels also very close to the Fermi en- 
ergy can be probed. Cs is evaporated from a thoroughly 
outgassed commercial getter source. The effect on the 
lifetime measurements by dosing a metal surface with a 
small amount of Cs has been thoroughly examined in 
the past^^. Using cross-polarized pulses no differences in 
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FIG. 2: 2PPE spectra of molybdenum recorded at varying 
photon energy; solid line, peak shift as expected for an unoc- 
cupied state; dotted line, peak shift as expected of an occupied 
state. 



the lifetime of volume states could be observed between a 
clean and a cesiated metal surface by means of TR-2PPE. 

Static 2PPE spectra of an uncesiated molybdenum sur- 
face recorded at different photon energies are displayed 
in Fig. 13 In accordance with conventional photoemission 
experiments, the position of the low-energy cutoff of the 
spectra is governed by the work function of the surface. 
The sharp high-energy cutoff on the other hand maps 
the Fermi edge of the metallic sample. 2PPE spectra 
from polycrystalline samples are typically characterized 
by a rather structureless spectrum due to the restricted 
spectral width accessible with the used excitation wave- 
length in the optical regime and due to a blurring of 
spectral features by the involvement of an intermediate 
state in the photoemission process. The appearance of a 
distinct peak (indicated by A) at -Exin ~ 0.6 eV (corre- 
sponding to an intermediate state energy oi E — Ep = 2.2 
eV) in the present example is insofar somewhat unusual. 
Rhodium, for instance, exhibits a completely structure- 
less 2PPE spectrum. According to the excitation pro- 
cess, the shift in the peak position by AE — hv as the 
excitation wavelength is tuned, indicates that the peak 
is due to an feature in the unoccupied DOS of Mo. In 
contrast, an occupied state feature would be character- 
ized by an energy shift proportional to A_E = 2hv (see 
dotted line in Fig.[21and the shift of the Fermi edge). Our 



measurements are performed on a polycrystalline sample 
and, furthermore, we do not observe any dependence in 
the appearance of the peak as a function of the excita- 
tion light polarization. We, therefore assign, this peak 
to an increase in the volume density of states in the cor- 
responding intermediate energy region rather than to a 
surface state. Our band structure calculations presented 



in Sec. IIV Al show indeed a drastic increase in the DOS 
at about 2.25 eV above the Fermi level (see Fig. 4). 

III. THEORY 

Within the framework of many-body perturbation 
theoryj^Si^ the inverse lifetime (decay rate) of an excited 
electron (hole) in a single-particle Bloch state il}qn with 
eigenvalue eq„ is given bjsi^ 



V = 2|ImAeq„|, 



(1) 



where Aeq„ is the many-body self-energy corrections to 
In the Dyson scheme, these corrections can be ap- 
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where Zqn is the renormalization factor related to the 
partial derivative of ASq„(a;) with respect to energy 
at cj = eq„. Here AEq„(ct;) is the matrix elements 
(i/'qn|AE(a;)|?/;q„) of the difference AE(a;) — S(a;)— w£^^ 
between the self-energy S(a;) and the exchange and cor- 
relation potential obtained in the local-density approxi- 
mation (LDA). Thus, the main problem in the lifetime 
evaluation is the calculation of the matrix elements of the 
self-energy i;q„(w). 

A widely used approximation to calculate the self- 
energy is the GWAi^ This approximation is formally 
the first-order term in the Hedin expansionS^ of the self- 
energy in power of the screened Coulomb interaction. In 
the simplest version of the GWA we consider here, the 
self-energy is defined by the product of the noninter- 
acting Green function G*' and the screened interaction 
W^ obtained within the random phase approximation 
(RPA)i22i224S Drawing attention only on momentum and 
energy dependence with the use of a shorthand notation 
q — (q, uj) in order to simplify the presentation, we can 
write the GW self-energy J^'^^ {cu) as^ 



^GW 



(9) 



(27r)4 



dkG"iq-k)W"{k). 



(3) 



Both G" and W'^ are built on the {ipun, Ckn} set. Physi- 
cally Eq. (PJ means that an electron, Cq, propagating in a 
system with the momentum g, creates a disturbance car- 
rying momentum k [see Fig.EJa)]. This fact is reflected 
in the formula by the interaction W°{k). Further, as de- 
picted by G°(g — k) in Eq. ©, the electron propagates 
with the momentum {q — k) right up to absorbtion of the 
disturbance. 



(_^)%L^ 




FIG. 3: A schematic representation of the GW self-energy (a) 
and the T-matrix contribution (b). The wiggly lines signify 
the RPA dynamically screened Coulomb interaction W^. The 
solid lines with arrows represent the Green function G". 



Using the spectral function representation of G" and 
W^° in Eq. Q , we can explicitly write down the imaginary 
part of the GW self-energy for, e.g., excited electrons 
{uj > ep) as 
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FIG. 4: Total and partial DOS in Mo. The total DOS is shown 
by solid line, the 4d DOS by dashed line, and the 5s+5p DOS 
by dotted line. 



where the T matrix is given by 



ImS^f (c.) = -E E E(^q«^k-q„'|i?k.) 



r(q) = W{q)P''{q)W{q) 



k n'>€F "ij 

X ImW°(k,w-ek-q„') 

X (Bkj|V'k_q„/?/jq„)e(w-ek-q„')- (4) 

Here tp is the Fermi energy, 9(x) is the step function, 
{-Bki(r)} is a set of appropriate Bloch basis functions, 
and W% are matrix elements of the RPA screened inter- 
action in this basis. In the present work, one-electron 
wave functions V'qm energies eq„, and i?ki are calculated 
by using the linear mufhn-tin orbital (LMTO) method.'^'* 
Having obtained the imaginary part of the self-energy, its 
energy-dependent real part can be found from the Hilbert 
transform. For details, we refer the reader to Refs. 1331 and 

El 

By making use of the T matrix that describes multi- 
ple scattering between an electron and a hole, we can go 
beyond the GWA in the calculation of the self-energy. In 
this case, we additionally take into account third- and 
higher-order terms in the Hedin expansion of the self- 
energy which constitute an infinite class of electron-hole 
ladder diagramsi^^^ This allows us to consider the fol- 
lowing scattering processes [see Fig. |2Ib)]. The e, elec- 
tron creates a disturbance which in its turn results in the 
creation of an electron(e)-hole(/i) pair. After that, the 
propagating electron repeatedly interacts (one or more 
times) with the hole during the e-/i-pair transport until 
this pair annihilation. The latter causes the disturbance 
which further is adsorbed by Cq at the end of the consid- 
ered process. 

Without gping into details which one can find in 
Refs. 23 and 24, we write the T-matrix contribution as 
an additional term to the GW self-energy of Eq. ||2J): 



1 + ^P°{q)W{q) 
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(5) 



The local electron- hole interaction W{q) is simply related 
with the first-order exchange correction to P" (q) account- 
ing for a single electron-hole scattering event. Note that 
Eq. © should be considered as a relation between the 
matrices 7ij(q, w), Wij{q,Lu), and P°(q, w). It is im- 
portant that the T-matrix contribution of Eq. (jSJ has a 
GW^-like formula. This means that in order to obtain the 
imaginary part of Eq„(w), we can use Eq. Q with Tlj 

instead of W". . 

''J 

The resulting self-energy in such a GW + T approach 
is given by the sum T,^^+^iu;) = 5]^^(w) + pLiuj). In 
calculations presented below, following Refs. l2335l and 
I12L we model the interaction W as the Fourier trans- 
form of the RPA screened interaction W'^ at zero en- 
ergy but with the 2/3 prcfactor obtained in Ref. |2J, i.e., 
Wij(k, cj) — |iy°(k, 0). Having obtained the self-energy, 
one can evaluate the lifetime with the help of Eqs. (Q) 
and ((D). In order to compare theoretical results with the 
experimental data (polycrystalline sample), we average 
the calculated lifetime over momentum for a given exci- 
tation energy. 



IV. RESULTS AND DISCUSSION 

A. Molybdenum 

Figure 01 shows the band structure of Mo and corre- 
sponding total and partial DOS. As mentioned before, 
the DOS has a wide (« 11 eV) d-band with a deep min- 
imum at Ep. Until ~ 5 eV above Ep the total DOS 
is almost completely determined by the unoccupied d- 
states and has a sharp peak at 2.25 eV. This peak can be 
identified with the experimentally observed unoccupied 
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FIG. 5: Momentum-averaged lifetimes r of excited electrons 
(e) and holes (h) in Mo calculated within GW and GW + T 
approximations. Inset, scaled momentum-averaged electron 
and hole lifetimes, t{E — -Ep)^, in Mo. 
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FIG. 6: Calculated and experimental momentum-averaged 
electron lifetimes r in Mo. Open circles are lifetimes obtained 
from a poly crystalline Mo thin film (thickness ~ 20 nm) , filled 
circles are lifetimes measured on a poly crystalline Mo sheet. 
Inset, scaled momentum-averaged electron and hole lifetimes, 
t{E- Epf, in Mo. 



peak A discussed in Sec. E] Below the Fermi level, the d 
states also dominate in the formation of the total DOS. 
At the same time, as it is evident from Fig. 0] the contri- 
bution of sp states is essentially larger for these energies 
than for energies above Ep. The total DOS shows two 
pronounced peaks (at -1.96 eV and at -3.4 eV) separated 
by a minimum at -2.52 eV. 

Figure \^ shows the calculated momentum-averaged 
lifetimes of electrons and holes. As follows from the fig- 
ure, the electron lifetime is shorter than the hole one. It 
is caused by both the considerable sp contribution below 
the Fermi level and the total DOS shape which differs 
from that above Ep (the available phase space for hole 
inelastic scattering is smaller). Nevertheless, as can be 
seen from the inset in Fig. |5l due to the fact that the 
total DOS has a relatively symmetric form with respect 
to Ep the electron and hole GW lifetimes demonstrate a 
quite similar behavior. Also, both lifetimes have a small 



knoll in the vicinity of the excitation energy of ~2 eV 
arising from the key features of the electron structure of 
Mo. As evident from the figure, the inclusion of the T 
matrix leads to a shortening of the lifetimes. An analysis 
of the ratio 7 = ImE<^^+^/Imi;<^^ has shown that in 
the vicinity of Ep (±0.5 eV) the T-matrix contribution 
is essentially less than the GW term. However, ImE'^ in- 
creases rather fast with increasing excitation energy and 
amounts on average to ~ 35% with respect to ImS''^ at 
|£;-£:f| «2eV. 

Addressing the role of d states in the formation of the 
energy dependence of the lifetime, it is worth emphasiz- 
ing that taking into account only d states in the evalua- 
tion of the e — h propagator P^ and the Green function in 
the GWA, one obtains similar results for both the elec- 
tron and hole lifetimes 1^ This means that the d-band 
plays a dominant role and only M states can be treated 
to quickly assess the imaginary part of the self-energy 
and the lifetime. 

In Fig. 1^1 the calculated momentum-averaged electron 
lifetimes are presented in comparison with the exper- 
imentally determined relaxation times of photoexcited 
electrons extracted from the cross-correlation traces. The 
highest excited state probed in the experiments was re- 
stricted by the photon energy to 3.0 eV. The transport of 
excited electrons out of the surface region into the bulk 
can act as an additional decay channel which reduces 
the measured lifetime of the excited states*^ In order to 
study the influence of transport effects on the experimen- 
tal results, we performed further experiments on the care- 
fully prepared Mo films suported by a dielectric MICA 
substrate. Different film thickness between 20 nm and 
100 nm were studied. As evident from Fig. El both, sheet 
data (solid circles) and thin film data (open circle), give 
identical results within the scatter of the data. The lack 
of any systematic deviation shows that the contribution 
of transport to the lifetime measured in the TR-2PPE 
experiment is negligible. 12*^ 

Regarding the comparison between experimental and 
theoretical data, inspecting Fig. |H| one can see that the 
T-matrix inclusion leads to better agreement between 
the theoretical and experimental results (especially for 
E — Ep > 1.0 eV). This becomes particularly evident in 
the inset of Fig. El which shows the deviation of experi- 
mental and calculated lifetime data from the Fermi-liquid 
model by Quinn and Ferrell"^^. Here, the center of gravity 
of the experimental scaled lifetime (w 25 fs eV^) is close 
to the GW+T curve for energies higher than ^ 0.7 eV. 
Below this value, towards Ep, the experimental lifetime 
goes down faster than the GW + T one. However, the ex- 
perimentally observed less rapid increase in the lifetime 
towards Ep in comparison to the Fermi-Liquid model in 
this energy region is qualitatively well reproduced by the 
theoretical results, which follow the slope of the data rea- 
sonably well. Some discrepancy for energies smaller than 
1.0 eV can be caused by both the approximations made 
for the self-energy (e.g., electron-phonon coupling can be 
important here) and secondary processes which effect the 




FIG. 7: Total and partial DOS in Rh. Total DOS is shown 
by solid line, the 4d DOS by dashed line, and the 5s+5p DOS 
by dotted line. 
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FIG. 8: Momentum-averaged lifetimes r of excited electrons 
(e) and holes (h) in Rh calculated within GW and GW + T 
approximations. Inset, scaled momentum-averaged electron 
and hole lifetimes, t{E — Ef)^, in Rh. 



experimental observed relaxation times i^ Note that this 
discrepancy can be eliminated by adding to ImS'^^+"'" a 
weak energy dependent term of '^ 5 meV. 



B. Rhodium 

The calculated band structure and DOS of Rh are 
shown in Fig. [7| Compared to Mo, the Rh d-band is 
more narrow {^ 7 eV). Almost all the d states are oc- 
cupied except for a small fraction which constitutes the 
sharp peak at 0.46 eV and pronounced edge of the d-band 
at 0.9 eV above Ep. At higher energies unoccupied bands 
show comparable contributions of sp and d states. Below 
E-p, the sp-states contribution to the bands is negligible 
in comparison with the d states. Another peak compa- 
rable with the previous one is situated at —2.58 eV. 

In Fig. |S1 the calculated momentum-averaged lifetime 
of electrons and holes in Rh is shown. On first glance, it 
may seem that the lifetimes of electrons and holes demon- 
strate a similar behavior and differ to only a small extent. 
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FIG. 9: Calculated and experimental momentum-averaged 
electron lifetimes r in poly crystalline Rh. Inset, scaled 
momentum-averaged electron lifetime in Rh. 



However, inspecting the scaled averaged lifetimes (see in- 
set in Fig. IHl , one can see that starting from an energy 
corresponding to the d-band edge, the energy dependence 
oi t{E — Ep)"^ for electrons can be fitted by a straight 
line with the slope of « 6.8 fs eV against « 2.2 fs eV for 
the hole scaled lifetime. Such a distinction is caused by 
the total DOS shape, the main part of which lies below 
Ep. Note that at energies below the d-band edge both 
the considered curves exhibit the scaled lifetime energy 
dependence similar to that in the case of Mo, but with 
essentially smaller values oi t{E — Ep)'^ due to the higher 
DOS.39 

As in Mo, the inclusion of the T matrix leads to a 
shortening of the lifetime. However, in this case, in the 
vicinity of Ep (±0.4 eV) the T matrix yields the largest 
contribution which amounts on average to ~ 40% in re- 
lation to ImE*^^. As the excitation energy increases, 
the contribution falls down. At \E — Ep\ « 3.0 eV, the 
ratio 7 has knolls running up to 1.13 for electrons and 
1.17 for holes. Calculations performed with taking into 
account only d states have corroborated^ the dominant 
role of the d-band in formation of the excited electrons 
dynamics in transition metals. 

In Fig. |S1 we show the experimental relaxation time 
in comparison with the calculated momentum-averaged 
electron lifetimes. Again, as it has been done for Mo, 
the highest excited state probed in the experiments was 
restricted by the photon energy to 3.0 eV. The lifetime 
measurements for polycrystalline thin films were not car- 
ried out in this case. Taking into consideration that the 
relaxation times of excited electrons in Rh are rather 
short over a large energy region, we have supposed that 
transport effects are negligible. Actually, the considered 
quantity remains below 7 fs down to an exciting energy 
E — EpK. 0.7 eV (see experimental points in Fig. |5J). To- 
wards the Fermi energy a steep increase in the relaxation 
time is obtained. This is the lowest onset energy in a life- 
time increase ever reported in TR-2PPE measurements 
of bulk states. 



Comparing the experimental and theoretical data 
shown in Fig. |^ one can distinguish the following main 
aspects. The inclusion of the T matrix ensures better 
agreement between the theoretical and experimental re- 
sults for E — Ep < 1.7 eV (see the inset), whereas for 
higher energies the experimental values are even above 
the GW curve. As in Mo, some discrepancy in the low 
energy region can be eliminated by adding to ImS'^^^-'" 
a weak energy dependent term, but the latter is equal to 
'--^ 10 meV in the case of Rh. 



C. Comparison of lifetimes in Mo and Rh 

In Fig. EH we compare the time-resolved 2PPE data 
as well as the theoretical results for the Ad transition 
metals rhodium and molybdenum. With a glance at 
the figure, practically over all the considered energy re- 
gion the electron lifetimes in Rh are shorter than those 
in Mo, especially below 2.0 eV. Presented as the ratio 

aW+T iGW+T 



'Mo 



/-, 



Rh 



(see inset in Fig. I10|l . the calculated 
GW + T lifetimes, which are consistent with the experi- 
mental data, clearly reflect this fact. Moreover, one can 
see that at energies above 1.0 eV, where agreement be- 
tween the experimental and theoretical results is rather 
good, this ratio falls down from 6 to unity (at ~ 3.5 eV). 
Further, the ratio becomes smaller than 1. By this com- 
parison we demonstrate that upon moving across the Ad 
series from Mo to Rh the changes in filling, width, and 
shape of the d-band have a strong effect on the electron 
lifetime. 

In principle, such a shortening of the lifetime is an ex- 
pected trend which can be qualitatively understood in 
a scattering theory approach (STA)iiii24^ Within this 
approach, the lifetime is considered in terms of energy- 
dependent matrix element |M(a;)p and the total DOS 
convolution. The former contains effects of screening, 
and the latter reflects peculiarities of one-electron en- 
ergy dispersion. As was shown in Ref. [la, for Mo and Rh 
a good accordance between the STA and GWA calcula- 
tions can be achieved with the energy-independent ma- 
trix element |M p. Furthermore, for excited electrons the 
matrix elements for Mo and Rh are similar. Therefore, 
an analysis of differences between the lifetimes in these 
transition metals can be reduced to an examination of 
features of the DOS. In addition, due to the dominant 
role of d states, only the partial d-DOS is of interest. For 
example, taking into account the density of d states av- 
eraged over the d-band width, as we move from Mo to 
Rh, one can obtain the electron lifetime decreasing by a 
factor of 3. Further, considering molybdenum as a repre- 
sentative of bcc metals from the center of the d-transition 
series and rhodium as a representative of fee metals from 
the end of these series, one can expect approximately the 
same shortening of the electron lifetime upon moving, 
e.g., from Nb to Pd (the Ad series) and from Ta or W to 
Pt (the bd series).^ 
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FIG. 10; Calculated averaged electron lifetimes r and exper- 
imental relaxation time in poly crystalline Rh and Mo. Inset, 
ratio of calculated electron lifetimes in Mo and Rh. 



V. CONCLUSIONS 

In conclusion, we have analyzed the strong influence of 
the transition metal Ad-hand structure on the relaxation 
dynamics of optically excited electrons. Comparing the 
obtained theoretical results on Mo and Rh with the ex- 
perimental data, we have shown that the T-niatrix contri- 
bution to the electron self-energy improves the GWA val- 
ues. Thus, we have found that taking into account mul- 
tiple scattering between an excited electron and a hole 
is important to achieve reasonable agreement with the 
TR-2PPE spectroscopy measurements. Next to the good 
agreement between experimental and theoretical lifetime 
data we were also able to assign a feature in the static 
molybdenum 2PPE spectrum to an enhanced electronic 
density of the excited state distribution in the calculated 
bulk band structure. Overall we believe that the present 
study is a good example of how the combination of so- 
phisticated state-of-the-art experimental and theoretical 
approaches enable a microscopic understanding of the ul- 
trafast femtosecond dynamics of electronic excitation in 
complex condensed matter systems. 
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